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Propagation attenuation rates for frequencies below 1 kc/s in the ELF region (30 c/s to 
3,000 c/s) were calculated from the spectra of atmospherics observed at widely-spaced sta- 
tions. Data are presented for east to west propagation under sunset approaching the eastern 
station. Under these conditions, the attenuation rates are about 1 db/1,000 km at 75 c/s and 
increase with increasing frequency, attaining about 3 db/1,000 km at 200 c/s. The attenuation 
rates observed seem to be consistent with a two-layered ionosphere model with its lower 
region 90 km above the earth. 



1. Introduction 

It is now accepted that radio waves in the ELK 
(extremely low frequency) region, extending from 
c/s to 3,000 c/s, propagate mainly in the waveguide 
bounded by the earth and ionosphere in the principal 
or zero-order mode. In recent years, a number of 
workers have investigated the propagation of ELF 
waves using atmospherics as signals, in lieu of more 
suitable man-made signals [Holzer and Deal 1956, 
Chapman and Macario 1956, Liebermann 1957, 
Pierce 1960, Smith 1960, and Watt I960]. 

In most of the work carried out to date, observa- 
tions were made at a single station of a large number 
of atmospherics having their origins in lightning 
discharges at different ranges from the receiver. The 
locations of individual lightning discharges were de- 
termined using direction finding systems or, in some 
cases, estimates of the storm locations were arrived 
at from meteorological information; thus, Chapman 
and Macario [1956] recorded waveforms at a single 
station while the locations of the lightning flashes 
were obtained from an atmospherics direction finding 
network. Propagation attenuation rates were then 
extracted from differences between the spectra of 
waveforms received at different ranges and from 
different sources. In this technique, some assump- 
tions regarding the constancy of the radiation spectra 
of the sources must be made. 

It would be of interest to determine propagation 
attenuation rates from waveforms of the same at- 
mospherics observed at different ranges; especially 
with the observing stations lying approximately 
along a great circle path passing through the origin 
of the lightning discharges. In this favorable cir- 
cumstance, fewer assumptions regarding the radia- 
tion spectra of the sources need be made. 

Such observations were made recently through the 
cooperation of the USAF who supplied atmospheric 
direction finding records from which the locations 



and times of occurrence of individual lightning dis- 
charges were determined. 

In this paper some preliminary results of atmos- 
pheric waveforms in the ELF band are present ed. 
Data are presented here resulting from fortunate 
alinements of storms and two recording stations 
approximately along great circle paths. The re- 
sultant attenuation rates are compared with theoreti- 
cal and other experimentally-determined values. 

The simultaneous recordings of atmospheric wave- 
forms of the vertical electric field were carried out at 
Boulder, Colorado; Fairbanks, Alaska; and Maui, 
Hawaii in the spring of I960. The equipment at 
each of the stations consisted of a 25 ft vertical 
antenna with an isolating cathode-follower con- 
nected directly to its base. The input impedance of 
the cathode-follower was approximately 400 meg. 
A coaxial transmission line connected the output of 
the cathode-follower with the amplifiers located in 
the recording building approximately 200 feet away 
from the antenna. The remainder of the equipment 
consisted of an amplifier, filters, oscilloscope and 
strip-film camera. The recording equipment and its 
performance were the saint 4 at each station. 

The overall amplitude versus frequency response 
of the recording apparatus, including a capacitor of 
100 fjifxi to simulate the antenna, is flat within 1 db 
between 30 c/s and 400 c/s and within 3 db between 
15 c/s and 1 kc/s (see fig. 1). The decrease in ampli- 
tude response at the higher frequencies, to eliminate 
interference from transmitters operating in the VLF 
region, was accomplished through the use of RC 
filters. The decrease in response at frequencies below 
45 c/s is a result of the time-constant of the interstage 
coupling capacitors and a reduction in voltage trans- 
ferred from the capacitive antenna to the cathode- 
follower. 

All of the receiving sites were located two miles or 
more from commercial power lines to minimize the 
reception of power line interference. Electrical 
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Figure 1. Amplitude versus frequency response of recording 
equipment. 

power was supplied to the recording equipment by 
diesel-electric generators. After taking these pre- 
cautions, it was found that the predominant inter- 
ference was received from the distant power lines. 
At Maui and Boulder, the 60 c/s interference (before 
filtering) amounted to about 5 mv/m and 10 mv/m, 
respectively. By comparison, the amplitude of 
received signals varied over a range of about 5 
mv/m to 100 mv/m. It was found advantageous to 
use a narrow-tuned filter to reject the 60 c/s power 
line interference. In plotting the spectra of the 
atmospherics, spectral components at 60 c/s were 
found by interpolating from components at fre- 
quencies either side of 60 c/s. 

Recordings were made simultaneously at each 
station on continuously-moving film, which was 
advanced past the cathode-ray tube at a speed of 
50 cm/sec. The film motion provided the sweep 
necessary to resolve the ELF components. 

2. Theoretical Considerations 

Let G(t,d) represent the amplitude variations 
of an atmospheric waveform in time, t, and at 
distance, d, such that 

| G(t,d) = for *<0 and t>r, (1) 

where r is the duration of the waveform. The 
Fourier spectrum of the waveform is then 

£(«,<©= C°G(t f d). e~ i(ai dt= ( r G(t } d)e' iut dt y 

Jo Jo 

(2) 

where co is the angular frequency. It is assumed 
that the source is equivalent to a vertical dipole and 
that the ELF waves of interest here propagate in 
the earth ionosphere waveguide in the zero-order 
mode. The spectrum amplitude of such a pulse, 
propagating in the zero-order mode, can be repre- 
sented by an expression of the form 



where A(u) represents the amplitude spectrum of 
the source function, a is the earth's radius in kilo- 
meters, and a(co) is the attenuation in nepers per 
unit distance. Equation (3) is valid at ELF pro- 
vided the sources are at least 1,000 km from the 
receiving antennas [Wait, I960]. It is also necessary, 
of course, that the observer is not near the antipode 
of the source. 

If the waveform of the same atmospheric is 
observed at two distances, d x and d 2l where d 2 ^>d iy 
it follows that 



exp [«(«) • (d 2 —d l )]-- 



\E(s*A)\ 



[sin dj/flTI* 
sin d 2 /aj 



(4) 



Using this expression, the attenuation rate as a 
function of frequency was obtained from the spectra 
of atmospherics observed simultaneously at distances 
d l and d 2 . Then the attenuation rate in decibels 
per 1,000 km of path length is given by 



a(ca)- 



f . , \E (w,di) | t n ! sin d 2 [al 1 3 

= 20 log : -r, / 7x 4 —10 log - — tt- ■ i r> 

L 6 | E(w,d 2 ) | & sin di/aj d 2 —d x 



(5) 
where the logarithms are to base 10. 

The amplitude spectra of atmospheric number 1 
observed at the three stations are shown in figure 2. 
The spectrum of highest amplitude represents 
observations made at Boulder, the shortest range; 
the spectrum of smallest amplitude represents 
observations made at Maui, the greatest range; while 



E(<*,d)mA(u)[a sin (d/a))- 1/2 e - 



x(co)d 



(3) 
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Figure 2. Amplitude spectra of atmospheric number 1 
observed at Boulder, Fairbanks, and Maui. 
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the intermediate values represent the Fairbanks 
observations. It is evident from these spectra that 
the frequency of peak amplitude decreases with 
increasing distance. For example, the frequency of 
maximum energy of atmospheric number 1 is about 
300 c/s at 4,500 km and 100 c/s at 9,600 km. It will 
be noted that the spectral width of the atmospheric 
narrows with increasing distance. The range of 
frequencies over which the attenuation can be cal- 
culated is limited to the width of the narrowest 
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Figure 3. Amplitude spectra of atmospheric number i. 
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Figure 4. Amplitude spectra of atmospheric number 



spectrum. The spectra and bandpass characteristic 
of the receiver indicate that the attenuation rates 
for the Boulder to Maui path should be valid over a 
frequency range of about 50 c/s to about 250 c/s. 

In a similar manner, the spectra of atmospherics 
numbers 2, 3, and 4, observed at Boulder, Fairbanks, 
and Maui, are given in figures 3, 4, and 5. 

In figure 6, the propagation attenuation rates in 
db/1,000 km derived from atmospherics numbered 
1, 2, 3, and 4 are plotted against frequency for the 
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Figure 5. Amplitude spectra of atmospheric number .',. 
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Figure 6. Propagation attenuation rates versus frequency 
derived from atmospherics 1-4 compared with calculated 
rates (see text). 
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Boulder to Maui pair. Attenuation rates derived 
from numbers 1,2, and 3 were recorded when sunset 
was at the eastern end of the path (1835 MST) and 
have minimum values of about 1.25 db/1,000 km 
near 75 c/s. Somewhat lower rates are observed at 
75 c/s for atmospheric number 4 which was recorded 
an hour later. It is tempting to attribute the lower 
attenuat'on rate observed at 1935 to the increased 
amount of darkness along; the path — although the 
occurrence of sunset along the path complicates the 
situation. 

3. Presentation and Discussion of Results 

During two, 2-minute recording intervals beginning 
1835 and 1935 MST on March 8, 1960, waveforms 
of five atmospherics were recorded at Boulder, 
Maui, and Fairbanks from storms in the Atlantic 
Ocean at approximately 39° N and 58° W. This 
information is tabulated below. 



Atmospheric 


Approximate range 
from source to 


Time of reception 

of atmospheric at 

Boulder 


No. 


Location 


Boulder 


Maui 




1 
2 
3 
4 


39°N, 59°W 
39°N, 65°W 
37°N, 61°W 
39°N, 58°W 


4500 
3500 
3500 
4500 


9600 
8700 
8800 
9600 


1835:52 MST 
1836:25 
1836:07 
1935:58 



The propagation paths from the sources to Boulder 
were in darkness. At the time of the first observation 
(1835 MST), ground sunset was approximately 600 
km west of Boulder and approximately 1,900 km 
west of Boulder during the second observation. 
In the first case, approximately 90 percent of the 
path from Boulder to Maui was in daylight; in the 
second case, approximately 75 percent of the path 
from Boulder to Maui was in daylight. 

The paths of propagation from the lightning dis- 
charges to Boulder and Maui, respectively, are 
approximately east-to-west and transverse to the 
earth's magnetic field. Waveform observations were 
also made at Fairbanks, Alaska and attenuation 
rates were calculated from the Boulder and Fairbanks 
data also, although this is a less favorable situation 
since the paths of propagation are not along great 
circle paths. 

As a matter of interest, calculated attenuation 
curves are also shown in figure 6. A two layer 
ionospheric model [Wait 1958] is employed; the 
height of the bottom of the reflecting layer is denoted 
h, and the effective conductivity parameter is 
w' r . At height d above this, the effective con- 
ductivity is assumed to jump to a new value of co'/. 
The curves shown in figure 6 correspond to the two 
cases where 



(i) co' r = 10 5 , h=90 km, d= co (upper calculated 
curve") 



(ii) a/ r = 10 5 , h=90 km, d~- 
curve) 

co'/ = 10 6 



= 20 km (lower calculated 



The first of these cases corresponds, of course, to a 
homogeneous ionosphere. The resulting attenuation 
should be proportional to (frequency) *. This would 
lead to a convex curve which is clearly not in accord 
with the data, The second of these curves has the 
required form; the ionosphere is thus behaving as a 
stratified conductor. For the lower frequencies, 
the upper well-conducting layer is becoming effective 
and the attenuation is reduced. Of course, other 
combinations of d and co" could be chosen to fit 
the facts, but it would always be necessary to have 
the effective conductivity increase with height in 
some manner. 

An average of the attenuation rates observed for 
the Boulder to Maui path is shown in figure 7 (solid 
line). There appears to be minimum in the atten- 
uation rate of about 1.0 db/1,000 km of path length 
near 70 c/s. The rate increases with increasing fre- 
quency attaining values of about 1.5 db/1,000 km 
at 135 c/s, 2 db at 160 c/s, etc, The increase in 
observed attenuation at frequencies below 70 c/s 
is doubtful; it is thought that the data are influ- 
enced by noise at frequencies below about 70 c/s. 
However, it is possible to construct a layer model 
which does have an attenuation increase as the 
frequency becomes very much lower. Such a model 
might consist of a three-layer ionosphere with para- 
meters co/, a/', and co/". If co/" is taken to be 
somewhat less than co/' and co/, a lower frequency 
will be reached where the attenuation is controlled by 
the outermost layer of conductivity co/" [Wait 
1957]. It is also possible that heavy ions may 
influence the results in frequencies below 100 c/s. 

Attenuation rates derived from the Boulder to 
Fairbanks and Fairbanks to Maui data are also 
plotted in figure 7. These data arc subject to ad- 
ditional limitations in interpretations; however, the 
rates are about the same as those presented for the 
largely sea-water Boulder to Maui path. In view 
of the theory for ELF propagation [Wait, 19601, the 
effects due to ground conductivity would be 
expected to be small if present. 

The propagation attenuation rates in figure 7 are 
compared with more complete data for day and for 
night by Chapman and Macario in figure 8. The 
present data plots between the daytime and night- 
time values reported by Chapman and Macario. 
Perhaps this is to be expected since the paths were 
partially dark at the east end, and in daylight to the 
west. It is of interest to recall that the Chapman 
and Macario data were recorded at ranges up to 
4, 000 km while the present data results from obser- 
vations at a minimum range of about 4,000 km and 
a maximum range of about 10,000 km. Although 
the present data are not very comprehensive, they do 
confirm attenuation rates resulting from other ex- 
periments using different techniques. In addition, 
the feasibility of deriving propagation data at ELF 
from spectral analyses of atmospherics is established. 
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Figure 7. Average of propagation attenuation rates versus 
frequenci/ for atmospherics 1-4 . 
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